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ABSTRACT 

The rising luminosity of the recent, nearby supernova 201 lfe shows a quadratic dependence with time during 
the first ss 0.5 — 4 days. In addition, the composite lightcurves formed from stacking together many Type la 
supernovae (SNe la) show a similar power-law index of 1 .8 ± 0.2 with time. I explore what range of power-law 
rises are possible due to the presence of radioactive material near the surface of the exploding white dwarf 
(WD). I summarize what constraints such a model places on the structure of the progenitor and the distribution 
and velocity of ejecta. My main conclusion is that the rise of SN 201 lfe requires a mass fraction X 5(> «3x 10~ 2 
of 56 Ni (or some other heating source like 48 Cr) distributed between a depth of ps 4 x 10~ 3 - 0. 1M S below the 
WD's surface. Radioactive elements this shallow are not found in simulations of a single C/O detonation. 
Scenarios that may produce this material include helium-shell burning during a double-detonation ignition, a 
gravitationally confined detonation, and a subset of deflagration to detonation transition models. In general, the 
power-law rise can differ from quadratic depending on the details of the event, so comparisons of this work with 
observed bolometric rises of SNe la would place strong constraints on the distribution of shallow radioactive 
material, providing important clues for identifying the elusive progenitors of SNe la. 
Subject headings: hydrodynamics — shock waves — supernovae: general — white dwarfs 



1. INTRODUCTION 

The use of Type la supernovae (SNe la) as cosmologica l 
distance indicators (Ries s et alJI 19981 IPerlmutter et al.lll999t) 
has brought attention to the theoretical uncertainties that re- 
main about these events. It is generally accepted that these 
SNe result fr om unstable thermonu clear ignition of degen- 
erate matter dHovle & Fowler! [i960) in a C/O white dwarf 
(WD), but frustratingly the specific progenitor systems have 
not yet been identified. The main three candidates are (1) 
stable accretion from a non-degenerate binary companion un - 
til the Chandrasekhar limit is r eached (Whela n & IberJl973h . 
(2) the mergin g of two C/O WDs dlben & Tutukovl Il984t 
Webbink 1984), or (3) accreting and detonating a helium 
shell on a C/O WD that leads to a prompt detonation o f the 
core (IWoosley & Weaver! [l994at iLivne & Arnettl fl995l) . A 
variation on the latter case is ignition triggered by a detona- 
tion i n an accretion stream (Guillo chon et alJ[2010tlDan et alJ 
2012). In addition, it is not known wheth er the incineration it- 
self p roceeds as a sub-sonic deflagration dNomoto et al.ll 197^ 
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1984) or deflagration-detonation transition (DDT; Khokhlov 



IWoosley & Weaver! Il994bl) . Single detonations of a 
sub-Chandrasekhar WD have been s hown to reproduce many 
features of SNe la dSim et al.1 120101) . but it is not clear how 
to ignite these cores without first detonating a helium shell. 
Each combination of situations has different implications for 
the velocity and density structure of the exploding WD, as 
well as the distribution of ashes. 

A potentially powerful method for constraining between 
these models are comparisons with the early-time behavior 
of SNe la, since this is when the surface layers of the WD are 
being probed by the observed emission. The recent, nearby 
SN 201 lfe is especially useful in this regard, since it was de- 
tected just fa 1 1 hours post expl osion when the lum inosity was 
merely ~ 10" 3 of that at peak dNugent et al.ll201 lh . Further- 
more, upper li mits on the luminos ity were placed ps 4 hours 
post explosion (Bloo m et al.ll2012l) . An interesting feature of 
the rise is a t 2 dependence for the luminosity up to pa 5 days 



post explosion. Furthermore, t here was no sign of the cooling 
of shock heated surface layers (Piro et al. 2010; Rabinak et al. 
1201 ll) nor interaction with a companion (Kasen 2010). This 
puts tight constrain ts on the progenitor radius of < 0.027? Q 
(Bloo m et al.ll2012h . demonstrating that it was very compact 
and consistent with a WD. Other individual SNe la have not 
been studied in the detail of SN 201 lfe, but the composite 
lightcurves formed from stac king many superno vae show a 
power-law index of 1.8 ±0.2 dConlev et al.ll2006l) . Although 
this is roughly quadratic, it could instead indicate some diver- 
sity in the rise. 

A f 2 dependence is consistent with a model in which the 
effective temperature remains fixed while the rad ius increases 
with time at constant velocity (Rie ss et alJI 1 999h . This seems 
unlikely to hold for a real supernova, since the effective tem- 
perature and photospheric radius can potentially change as 
the ejecta expands and its density d r ops. In the Supple- 
mentary Information of Nuge nt et alJ (1201 lh . a sing le-zone 
mode l is described (using arguments similar to that in I Arnettl 
1982), estimating the luminosity from expanding ejecta that 
is heated by 56 Ni decay. For a explosion energy E SD and to- 
tal stellar mass M*, one can define a characteristic velocity 
of V = {2E sil /M*) l l 2 . Using a constant opacity n, this gives a 
luminosity of 



L(t): 



2ir cVX 56 e 56 t 2 



(1) 



where c is the speed of light, is the 56 Ni mass fraction, 
and eg6 = 3.9 x 10 10 erg g" 1 s" 1 is the radioactive heating rate 
per unit mass. This result demonstrates that radioactive heat- 
ing can in principle also provide a f 2 dependence. What this 
single-zone model cannot answer is what is the required depth 
of the radioactive material. In addition, the power-law depen- 
dence on time may change depending on gradients in density, 
velocity, and distribution of radioactive isotopes. 

In the following work I generalize this single-zone model 
to a one-dimensional calculation that includes the above men- 
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tioned complications. In ^2j I summarize the main equations 
used to describe the dynamics and thermodynamics of the 
heated, expanding ejecta. In ^3] I discuss the time-evolving 
energy density of the expanding layers, and show why heat- 
ing from radioactive decay dominates over the shock heat- 
ing that has been the focus of previous work dPiro et alJl2010b 
Rabinaket al. 2011). I also explore the depth of the thermal 
diffusion wave as the surface layers expand, which shows the 
depth in the WD that the early lightcurve is probing as a func- 
tion of time. In 5jU I calculate the time-dependent luminosity, 
and explore what range of power-law scalings are expected for 
the rise. In $5j I compare these calculations with the obser- 
vations of SN 201 lfe to constrain the distribution and amount 
of radioactive material. Finally, in ^6]l summarize my results 
and discuss what detailed modeling can do for the understand- 
ing of the outer ejecta of SNe la. 

2. DYNAMICS AND THERMODYNAMICS OF THE 
EXPANDING STAR 

For this study I concentrate on the plane-parallel surface 
layers of the exploding progenitor star. This does not mean the 
expansion itself is plane-parallel, but merely that all the ejecta 
originates from roughly the same radius. This simplification 
is an accurate representation of the outer material on a WD. 
In the Appendix, I summarize the scalings for a non-plane- 
parallel treatment. 

Variables of the pre-expanded progenitor star are denoted 
with the subscript 0. The profile is assumed to be a polytrope, 
with 



Po=K Pt 



1+1 /n 



(2) 



where in the case of non-relativistic electrons n = 3/2 and 

K = 9.91 x I0 l2 p e 5 ^, and for relativistic electrons n = 3 

and K= 1.23 X I0 15 p e 4 ^, where p e is the molecular weight 
per electron and K is in cgs units. For a constant flux, 
ideal gas dominated, plane-parallel atmosphere n = 3 and 

K = 6.1 x 10 13 gq^T^g 5, where g = I0 9 gg cm s" 2 is the sur- 
face gravity and T s s = 10 5 ?rff s K is th e effective tempera- 
ture of the photosphere dPiro et alJl2010h . Depending on the 
progenitor model, or heating from sub-sonic burning during 
the explosion, a hig her T e s may be a more realistic choice 
dRabinak et al.ll201lh . 

For a typical Chandrasekhar-mass progenitor, n = 3/2 is 
most relevant since at « 0.5 day post explosion (when 
SN 201 lfe was first observed), the diffusion wave is at a 
depth where the progenitor star equation of state is set by 
non-relativistic electrons. Nevertheless, I keep the poly- 
tropic index as a free variable since it can vary in other 
situations. For example, if the WD was ignited from a 
detonating heli um shell, the layer may have a convective 
profile instead (iBildsten et alJ 120071: IShen & Bildstenl 120091: 
iWooslev & Kasenll201 lUShen et al.11201 11 ). 

The velocity of the shock has a gradient with density, which 
rises toward the surface according to dMatzner & M cKee 
1999), 



Vs(po) 



1 0.79 ( — 



1/2 



(3) 



where /?„ is the WD ra dius and J3 = .19 for a radiation pres- 
sure dominated shock (Sakurai 1960). I keep /3 as a free pa- 
rameter, so that my results are general with respect to veloc- 
ity gradients that are from non-shock sources. Typical val- 



ues for equation (O are E SD w 10 51 erg and the mass and ra- 
dius of a WD near the Chandrasekhar limit, but it may also 
be worth considering a lower mass WD in light of the pure 
de tonation simulatio ns of sub-Chandrasek har explosion mod- 
els (ISim et al.ll2010l) . In lPiroetal] (120101) . we focused on the 
shock from a DDT, and wrote the velocity gradient as 



V s ( Po ) = V (po/p') 



-13 



(4) 



where V' and p' are set by where the detonation fails and 
a shock runs away, heating the surface of the star. This 
roughly gives p' f=s 2 x 10 6 g cm" 3 and V' « 8 x 10 8 cm s" 1 
dPiro et alJ l2010). For the present work I quote analytic re- 
sults in terms of p' and V' since it gives the cleanest solu- 
tions. But for the numerical factors, I substitute p' = M*/R\ 
and V = 0.79(E sn /M*) 1 / 2 , since using E sn and makes it 
easier to compare with other theoretical work and observa- 
tions. 

In this framework, one can think of the surface layers as 
a series of shells, each labeled with an initial density po and 
moving with a final velocity of V(po) « 2V s (Po) (this factor 
of 2 is the effect of pres sure gradients causing acceleration, 
Matzner & McKee 1999). For a plane-parallel configuration, 
mass conservation as the shell expands gives 



p(p 0l t) = po 



R* 



V(p )t 

1/m 



Ho 



AV(p )t 



(5) 



where Hq = Po /pog = Kp jg is the thickness of the layer, 
which I set to be the pressure scale height, and the velocity 
gradient of the layer is 



opo on) 



l + l/n 



V(po). 



(6) 



This expression is found by using the equation of hydrostatic 
balance, dPo/dro = -pog, and equation ©. Thus AV(po) is 
smaller than V{po) by a constant factor of (3/(1 + 1/n) w 0.1 1 
(for n = 3/2 and /3 = 0.19). 

The thermal evolution of the expanding layer is set by the 
first law of thermodynamics, 



Tds = <l I — ) +Pd (-) = -dU+-Ud (- 

P ) \PJ P 3 \p 



(7) 



where s is the specific entropy, U is the energy density, and 
the right-hand side comes from assuming a radiation domi- 
nated energy density, so that P = U/3. Changes in entropy 
come from nuclear heating and radiative losses, so this can be 
rewritten as 



1 dU 4 d ( 1 
~p~dt + l> U dt \p 



(8) 



where the partial derivatives in time are evaluated at a fixed 
mass shell, labeled by po- I assume that the 56 Ni can poten- 
tially vary with depth, which is modeled with a characteris- 
tic density p^ and a power-law index \- Such a power-law 
choice is not physically motivated, but is made simply to al- 
low the deposition to vary, while still resulting in self-similar 
solutions. I ignore changes to due to decay, since I am 
focusing on times earlier than the 56 Ni half-life of 6.077 days. 

It is also possible that other radioactive isotopes could 
be present and be powering the early-time lightcurve. In 
particular, 48 Cr has a 21.56 hr half-life with nearly 100% 
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electron captures to the 1 + excited state of 48 V, which is 
followed by a cascade that emits an energy of 0.42 MeV. 
The 48 V has a 15.973 day half-life, and an effective en- 
ergy deposition of 2.874 MeV. Therefore the first step de- 
posits w 7.5 x 10 1() erg g" 1 s" 1 , and the second step deposits 
w 2.9 x 10 10 erg g" 1 s" 1 . These are not too different than eg6, 
but may introduce some qualitative changes. The radioac- 
tive isotope 52 Fe may also be present, but its short half-life 
of 8.275 hr means that it is negligible at all but the earliest 
times. For the main analysis of the paper, I concentrate on 
56 Ni, but one should keep in mind that the factor X^ese can 
stand for a more complicated mixture of radioactive isotopes. 

The derivative dL/dM r is the radiative cooling rate per unit 
mass, where 



L = - 



4irr 2 c dU 
3np dr 



(9) 



Throughout my main analysis I assume that k is independent 
of density and temperature, such as for an electron scattering 
opacity. This is roughly accurate for the hot surface layers, 
but in the Appendix I summarize the self-similar scalings that 
result from a more general opacity law. 

3. ENERGY DENSITY AND DIFFUSION DEPTH 

Using the framework summarized in the previous section, 
one can in principle solve equation (H)) for U(po,t). This 
provides the entire dynamic and thermodynamic evolution 
of each mass shell (labeled by po) as a function of time. 
For the present work I am focused on the early lightcurve 
(« 0.5-4 days post explosion) in the limit where heating from 
56 Ni decay dominates. The strategy is then to assume that the 
nuclear heating is much greater than radiative cooling, so that 
dhj dM r can be ignored in equation ([8]). Once this simplifica- 
tion is made, the equation can be easily integrated to find 



U(po,t)-- 



l + l/nRlKp' M /"e 56 X 56 



16(3 



gV' 3 t 2 



£o_ 

P56 



l+l/n+3/3 



(10) 



where there is a factor of ss 1 /8 included from taking V « 2V S . 
This solution implicitly assumes that the energy density de- 
posited by the passing sho ck wave, which w as t he main source 
of hea ting considered by iPiro et ail (120101) and Rabina k et alJ 
is negligible. To show that this is a reasonable approx- 
imation, consider a shock traversing a density po- The energy 
density deposited by that shock is 



Ush.O 



■ 1 



PoVf 



(11) 



This energy density then adiabatically decreases as the layer 
expands, giving 



U sh ( Po ,t) = U s 



sh.O 



p(Po,t) 



psh 



(12) 



where 7 = 4/3 is the adiabatic exponent for radiation- 
dominated material and p^ = po(7 + l)/(7~ 1) = Ipo is the 
compressed, shocked density. Since t/ sh oc r 4 (from eq. fTOl ) 
while U oc f" 2 (from eq. [ (fTOl ). the shock energy density 
falls off much faster and is quickly negligible in comparison 
to the radioactive heating. I compare the two energy densi- 
ties in Figure 1 at three different snap-shots in time. This 




10 5 10 6 
Po (g cm- 3 ) 

FIG. 1 . — Comparison of the energy density due to radioactive decay given 
by equation (Tp) (dashed lines) for X 56 = lfT 3 , 1(T 2 , lfT 1 , and 1.0 (from 
bottom to top in e ach panel) with the energy density from shock passage 
given by equation (T2J (thick, solid lines). The supernova parameters are 
£ sn = 10 51 erg, M* = IAMq, andff* = 3 X 10 s cm. The vertical dotted lines 
show the depth of the thermal diffusion wave at each time (eq. 1151 ). Ra- 
dioactive decay dominates the energy budget for > 10 at 2.4 hrs after 
shock breakout and for > 10~ 2 at 7.2 hrs. 



shows that the radioactive heating dominates at m 1 day un- 
less X 56 < 10" 3 . 

Also plotted in Figure 1 is the depth of the thermal diffusion 
wave as the surface layers expand and cool (vertical dotted 
lines). The timescale for a given layer to cool is 



fdiff ; 



3kp 



(Ar) 2 



(13) 



By substituting equation (O into equation ( fT3l and setting 
fdiff = t, I solve for the diffusion depth as a function of time, 



P0,diff(0 



2{\ + \/n) V'cg 



1 i/(i+i/n+/?) 



3/3 nRlKp' l+x l" 



,2/(1+1 /n+/3) 



(14) 



For n = 3/2, j3 = 0. 19, and using the equation of state for non- 
relativistic electrons, this gives 

po,diff(f) = 2.5 x 10 k 02 t 51 M lA K g5 f day gem , 

(15) 

where K0.2 = k/0.2 cm 2 g" 1 , £51 = E sn /1Q 51 erg, Mi 4 = 
M»/1.4M , R S5 = #*/3 x 10 8 cm, and f day = t/l day. This 
density is roughly consistent with what one would expect for 
material dominated by non-relativistic electrons. To empha- 
size the depth in the WD this corresponds to, I estimate the 
diffusion mass as Mdiff ~ (2/5) x 47r/? 2 po.diff#o(po.diff) (where 
the factor of 2/5 is appropriate for an n = 3/2 polytrope), so 
that 



MdiffCO 



1.5 X 1U K 02 h 51 M lA K S5 t dRy MQ. 



(16) 
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FIG. 2. — The solid line marks the position o f the thermal diffusion wave 
as a the surface layers expand using equation 1 1 6i . This shows the depth 
in the WD from which the luminosity is originating as a function of time. 
Along the solid line I have plo tted filled circles at the time of the PTF ob- 
servations I Nugent et al. 2011) an d an o pen circle showing the time of an 
upper limit constraint (Bloom et al. 2012). The horizontal dashed lines show 
the mimimum helium shell mass necess ary for detonatio n, which are each 
labeled by their corresponding core mass i Fink et al. 2010). For example, for 
a 1.125Mq C/O core a ().039Mq helium shell will successfully detonate, as 
will any larger mass helium shell. 

The most striking feature of these results is that 56 Ni must 
be present « 1O~ 2 M below the surface at « 1 day post- 
explosion. The earliest detection of SN 201 lfe is at « 11 hr 
(assuming a f 2 rise), which would require 56 Ni at a depth of 
merely w4x 10" 3 M Q . 

In Figure |2j I plot the depth of the diffusion wave be- 
low the WD surface as a function of time using equation 
( TT6I ). The filled circles in dicate the time of the PTF observa- 
tions dNugent et alj|201 ll) . Since the time of explosion of SN 
201 lfe is not exactly known, I assume it is 55796.687 MJD , 
which was derived assuming a f 2 rise (iBloom etal.ll2012h . 
Also plotte d as an open circle is the time of the upper limit 
reported by IBloom et ail (120121) . 

As a comparison, I also show the minimum helium shell 
masses needed for a surface detonatio n (dashed lines, from 
iShen & Bildsteri2009tlFTnk et al.l2010l) . which are labeled by 
their corresponding core masses (a smaller C/O core requires 
a larger helium shell mass for ignition). Such surface helium 
burning can ignite the core and create a double-detonation ig- 
nition. This comparison between the shell masses and the 
depth of the diffusion wave in SN 201 lfe indicates that a 
double-detonation ignition may provide radioactive ashes at 
the correct depth needed for powering the early-time rise. 
I further discuss the double-detonation scenario, along with 
other SNe la progenitor models, in $6] 

4. RISING LUMINOSITY 

I next estimate the luminosity expected from a 
radioactively-dominated rising supernova, which will 
allow constraints to be placed on the amount and gradient of 



56 N near the WD surface. Combining equations (O and JTOb . 
the luminosity from a given shell is 



Upo,t) = 



1 + 1 jn + 3/3+ x 4ircV 'e 56 X 56 t 7 



3 k 



po_ 

P56 



(17) 



From this result, a number of important features can already 
be identified. In the limit that \ an d P are small, all dependen- 
cies on pq are negligible, so that L ~ cV'^es^/n, consistent 
with the single-zone result given by equation (Q}. 

To derive the luminosity as a function of time I substitute 
Po diffW from equation ( fl4l > in for p a in equation ( fT7] i, 



1(0= 



l + l/n + 3/3+x 



P 
3 k 



2(l + l/«) gV'c 



3/3 KRlKp' M / n 
r 



x-0 

]-] ll-;S 



AttcV' 656^56 ( p' ^ X f 2(l+l/B+ x )/(l+l/,i+/3) 
P56 



(18) 



Again, when (3 s=s \ ~ 0, the power law simplifies to L oc t 2 . As 
the gradient in V s becomes larger and /3 increases, the power 
law decreases below 2. Conversely, a larger deposition index 
X can increase the power law, since this means that more 56 Ni 
is being probed as the diffusion wave moves deeper. 
For the case of n = 3/2, j3 = 0.19, and \ = 0, the luminosity 



L(0 = 2.1xl0 42 X 56 ^° 2 90 £ - 45 



0.10 1.8 -1 
8.5 May er S S ■ 



(19) 

For SN 201 lfe, a lum inosity of » 10 40 erg s" 1 is seen at 
w 1 1 hr (Nugent et al. 201 1). From the estimates given here, 
this requires a mass fraction of X^e, ~2x 10~ 2 (although in 
the next section I make a more quantitative comparison to the 
observations). Such an amount of 56 Ni is consistent with the 
energy density of the plasma being dominated by radioactive 
decay as discussed in Figure 1 . 

For my fiducial model given by equation ( fl9l ), the rising lu- 
minosity scales like f 18 and not t 1 . This is inc onsistent with 
the p ower-law index of 2.01 ± 0.01 found by dNugent et alj 
1201 1 | ). but within the measurements of a broader range of SNe 
la dConlev et alj|2006l) . This then begs the question, is an in- 
dex of 2 universal across all SNe la, or can the power law 
vary around 2 depending on the gradients in velocity and ra- 
dioactive heating for a given particular event? The work here 
argues via equation ( fT8l that a range of power-laws are possi- 
ble. The fact that the best studied SNe la shows an index very 
close to 2 may be evidence that it is universal, but remember 
that the PTF data is in g-band. A detailed comparison with 
the bolometric rises of SNe la is necessary to understand how 
diverse the power-laws can actually be. 

If the bolometric luminosity does indeed follow a simple 
f 2 scaling, then the 56 Ni deposition must obey x 1=3 P- This 
cancels out the dependencies due to the position of the thermal 
diffusion wave (essentially the entire bracketed term in eq. 
fTSl is set to 1), giving the simple dependencies that can be 
derived from a single-zone analysis, 



L(f) = 2.7 x 10 42 X 56 



^^EfMlftlyKgs- 1 , (20) 

where now it is understood that Xs^.p^ is the mass fraction of 
56 Ni at po = 10 6 g cm" 3 , which corresponds to the depth of the 
diffusion wave at w 10 hrs (eq. fT5l ). 
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5. COMPARISONS TO SN 201 1FE 

The effective temperature of the supernova during these 
early times may be complicated due to an opacity that depends 
on a mix of heavy elements. With this caveat in mind, I esti- 
mate the effective temperature assuming a constant, electron- 
scattering opacity. The effective temperature is 



T ef f=T/r 



1/4 



(21) 



where T = (U /a) x l A is the local temperature and r » npAr is 
the optical depth, resulting in 



'eff : 



l_ + l/ne 56 X 56 y /4 fp YWpo\ 
4/3 anV) {p 5 J [p'J ' {ZZ) 



X/4 



/3/4 



Evaluating this at the diffusion depth provides an estimate for 
the time evolution of the effective temperature. For n = 3/2, 
/3 = 0.19, and x = 0, 



r eff (f)=i.6xioxr^ 



1/4-0.28 F -0.11 j^O.10 



56 ^0.2 



0.03.0.05 rv- 
'day ^ 



(23) 



This shows that the radioactive heating balances the expan- 
sion to give an effective temperature that is nearly constant 
with time. 

The photospheric radius r p h of the expanding ejecta can be 
estimated with relationship L = 47rr 2 h <7sB T^ff' where ctsb is the 
Stefan-Boltzmann constant. When the diffusion wave is at po, 
the photospheric radius is therefore 



r p h(Po,t)-- 

which for n = 3/2. 
tion ( TT5l l. results in 



l + l/« + 3^+xl 1/2 
3(l + l/n) 

P = 0.19, and X 



4V't ( £° 
p' 



(24) 



0, and using equa- 



r ph(f ) = 2.2 x lO'^^YW^W cm. 



' 1.4 



l 8.5 'day 



(25) 



Using this photospheric radius, and assuming black body 
emission with effective temperature given by equation ( 1231 1. 
I calculate the g-band absolute magnitude in Figure [3] Plot- 
ted as solid circle s are the observations summarized by 
iNugent et al.1 (1201 ll). along with an arrow indicating the up- 
per limit from Bloo m et al.1 (J2012). For this comparison a 
choice of explosion time is necessary, for which I use 55796.6 
MJD. This is not meant to maximize the fit between theory 
and data, but merely to provide a useful comparison. Al- 
though « 3 x 10~ 2 is favored, this model overpredicts the 
early lightcurve and underpredicts the late lightcurve, indi- 
cating that a gradient in 56 Ni deposition is needed to make a 
better match. This is not surprising since these models corre- 
spond to L oc f 1 8 and not f 2 . 

Motivated by this discrepancy between the observed 
lightcurve and that predicted for a constant distribution of 
56 Ni, I instead consider the case \ = ft = 0.19, which was 
found to give a f 2 bolometric luminosity (eq. ll20l ). The cor- 
responding effective temperature is 



T eS (t)= 1.6 x10 4 X^k^% 



-0.09,0.10 
1.5 



f d u a ; u K,(26) 

where again I emphasize that X^, P6 is the mass fraction of 
56 Ni at a depth of po = 10 6 g cm" 3 . The difference between 
equations (l23l and (l26b may not seem appreciable, but re- 
member that L oc T* s , so that small changes in exponents can 
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Fig. 3- 

isn = 10 51 
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Z56 = 3 X 10 3 to 1, as labeled on each curve. Although Z56 Ri 3 X 10 2 is 
favored, this model with a constant distribution of 56 Ni overpredicts the early 
lightcurve and underpredicts the late lightcurve. This indicates that a gradient 
in 56 Ni deposition is needed. 
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FIG. 4. — The same as Figuref3] but with a ,6 Ni gradient set to x = /9 = 0. 19 
instead. The model with -?56,p 6 ~ 3 X 10~ 2 provides a reasonable fit to the 
observed lightcurve. 



make an important difference. The photospheric radius is es- 
sentially unchanged from equation (l25l l, with only a small dif- 
ference in the prefactor. 

The resulting lightcurves for the x = 0.19 model is plotted 
in Figure |4] which demonstrates a better fit to the data. The 
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necessary gradient in 56 Ni implies that m 3 x 10~ 2 at a 
depth of Mdiff ~ 4 x 10~ 3 M Q and this varies up to X 5f> w 5 x 
10~ 2 at a depth of Mdiff ~ 0.2M Q . The total integrated mass 
of radioactive material needed to produce the lightcurve over 
» 4 days is roughly « 10~ 2 M Q . 

6. CONCLUSIONS AND DISCUSSION 

I investigated the impact of radioactive heating on the early- 
time rise of SN la lightcurves. A mass fraction X$(, «3x 10~ 2 
of 56 Ni at a depth w4 x 10" 3 -0.1M Q from the progenitor sur- 
face is needed to produce the luminosity seen from SN 201 lfe 
during its first « 4 days (although, as discussed in §2, other 
heating sources such as 48 Cr could just as well explain the 
rise). A model in which there is a velocity gradient set by 
the passage of a shock, but with a constant deposition of ra- 
dioactive material, gives a t 1,8 power law. This appears to be 
inconsistent with SN 201 lfe, but is not ruled out by the stack- 
ing of many nearby SNe la (IConlev et a l. 2006). The shape 
of the SN 201 lfe lightcurve is better fit when the gradients in 
velocity and radioactive material are similar, with \ ~ P, as 
is shown in Figure |4] In the future, a comparison between the 
bolometric rising luminosity of SN 201 lfe and this work us- 
ing equation cTTSb would provide important constraints on the 
distribution of radioactive material near the WD surface. 

A crucial question for future observations is whether the 
rise of SNe la obey a universal power-law or if they vary from 
event to event. My work shows that the rise should depend on 
the particular gradients in velocity, density, and deposition of 
radioactive elements, which, although it may be close to f 2 , 
should not necessarily always be the same. If the power law 
is indeed found to be universal, this would argue that different 
physics than what I am exploring here is determining the rise. 
For example, an opacity effect that I have not included could 
enforce a fixed T e fi and constant expansion velocity. 

The total integrated mass of radioactive material needed to 
produce the SN 201 lfe lightcurve over w 4 days is w 10~ 2 M Q . 
Single detonation models of Chandrasekhar and even sub- 
Chandrasekhar-mass WDs find distributions of ashes that are 
fairly well-stratified, and do not show 56 Ni or 48 Cr near the 
surface. This is perhaps not that limiting of a constraint; 
Chandrasekhar-mass single detonations are disfavored be- 
cause th ey cannot produce the observed intermediate mass el- 
ements (Filippenko 1997), and although sub-Chandrasekhar 
single detonations match the nucleosynthesis generally seen 



from SNe la (ISim et al J 120101) . it is not known how to ignite 
such an object without a helium shell. 

Conversely, the explosive ignition of a helium shell in 
the double-detonation scenario can produce shallow radioac- 
tive material. Indeed, Figure 2 shows that the depth of 
this material is roughly the same a s the minimum helium 
shell masses nee ded for detonation (IShen & Bildstenl 120091; 
iFink et all 120 10)) . The total amount of radioactive material 
needed is also qualitatively similar to the nucleosynthesis of 
such events (Fink et al. 2010). DDT models can produce 56 Ni 
near the WD surface, but potenti ally only in the mo st strongly 
mixed, off-center deflagrations (M aeda et alJl2010l) . In DDT 
models with many ignition points that have fairly stratified 
ashes, radioactive elements are not present near the surface. 
A gravitationally confined detonation (GCD) is in a sense 
just a more extreme, off-center version of the DDT models, 
and it too produces iron-peak elements near the surface when 
a bubble unstably rises and breaks at the top (Mea kin et al.1 
2009). Finally, a more speculative idea is that some radioac- 
tive elements are synthesized near the surface by ff -mode heat - 
ing during the pre- explosive convective phase (Piro 2011). 
iNugent et all (120111) report the presence of O, Mg, Si, S, Ca, 
and Fe in the spectra of SN 201 lfe at early times. Many 
of these elements are potential ashes from the scenarios de- 
scribed above, and a more detailed comparison may help dis- 
criminate between them. 

Even though this work argues for the presence of radioac- 
tive material near the surface of the WD, a potential problem 
is that if the abundance of iron-peak elements is too high, they 
tend to produce colors that are too red and sp ectra that are 
incon s istent with nor mal Type la supernovae (iKromer et al.l 
120 10t ISim et al.l 1201 lb . Although it should be noted that 
this difficulty may be partially alleviated for larger mass 
WDs that have smaller helium shells, and may also depend 
on careful consideration of the burning in the helium shell 
(IWoosley & Kasen1l201 ll) . Combining detailed modeling of 
the early lightcurve rise with comparisons to the peak colors 
and spectra provides competing limits on the mass fraction 
of iron-peak elements, and therefore together should result in 
tight constraints on the composition of the outer layers. 

I thank Lars Bildsten, Luc Dessart, Dan Kasen, Peter Gol- 
dreich, Christian Ott, and Ken Shen for helpful comments 
and discussions. This work was supported through NSF grant 
AST-0855535 and by the Sherman Fairchild Foundation. 



APPENDIX 



A. GENERAL SELF-SIMILAR RISING LUMINOSITY SOLUTIONS 

Here I consider a wider range of self-similar solutions for the rising luminosity. This highlights the expected changes in the 
time-dependence as various details are added to future models. 



A. 1 . General Opacity Law 

For the majority of this work I have assumed a constant opacity, consistent with electron scattering. To make better comparisons 
with observations will require complete calculation with detailed opacities. If helium is present near the WD surface (such as in 
the double-detonation scenario), it may recombine in expanded, cooled layers. Metals lines would provide a strong opacity in the 
UV. In light of this, I consider a more general opacity law 



K = Ko(fT b . 

The self-similar solution for such a case results in L oc t x (analogous to eq. ifTSl ). where 



A= (l + ia + b/l) 



l + l/n + x + xb/2 



l + l/n + P+(l + l/n + 3f3)(a + b/4) + X b/4 



(Al) 



(A2) 
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For example, for a pure C/O composition Rabinak & Wax manl d201 ll) use k = 0.66(r/10 5 K) 127 cm 2 g -1 . Combining such an 
opacity law with n = 3/2, f3 = 0.19, and \ = results in A = 1.71 for the power-law exponent. When instead x = P — 0.19, then 
A = 1.98, which shows that a roughly quadratic luminosity increase is fairly robust as long as the nickel deposition is increasing. 

A. 2. Non-Plane-Parallel Solution 

In the continuity relation, equation @, I made the plane-parallel assumption that essentially all the exploding material came 
from roughly the same radius of This is a reasonable approximation for a compact progenitor like a WD. Conversely, for a 
more extended progenitor continuity becomes 



p(po,t) = po 



A) 



V(po)t 



\AV(p )t 



(A3) 



where ro is the initial radius for a given shell of material. For a polytrope, this scales ro oc p at depths well below the surface 
of the star. Completing the self-similar analysis with a constant opacity, the diffusion wave now has a time-dependent position 
of 

Po.d.ffWocf^ 1 - 1 /^). (A4) 

This generally gives a larger exponent than the plane-parallel case in equation ( TBl i. showing that the diffusion wave traverses 
through the star more rapidly in this case. The time dependence of the luminosity is 

Locf 2(l-l/„ +x )/(l-l/„ + « (A5) 

For n = 3/2, j3 = 0. 19 and x = 0, this results in an exponent of 1 .27, which is much below 2. This means that as non-plane-parallel 
effects become more important, the luminosity will begin to flatten. Also note that this new exponent is much more sensitive 
to \, since setting x = P increases the exponent all the way up to 2. These non-plane-parallel solution may also have use in 
investigating Type Vole supernovae at early enough times that the li ghtcurve is still ris ing, but not so early that the emission if 
dominated by the shock-heating of the surface layers (as studied by Nakar & Sari 2010). 
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